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Abstract 

We study the influence of the velocity dependence of friction on the escape of a Brownian particle 
from the deep potential well (Eb 3> ksT, E^ is the barrier height, ks is the Boltzmann constant, 
T is the bath temperature). The bath-induced relaxation is treated within the Rayleigh model 
(a heavy particle of mass M in the bath of light particles of mass m <C M) up to the terms of 
the order of 0(A 4 ), A 2 = m/M <C 1. The term ~ 1 is equivalent to the Fokker-Planck dissipative 
operator, and the term ~ A 2 is responsible for the velocity dependence of friction. As expected, 
the correction to the Kramers escape rate in the overdamped limit is proportional to A 2 and is 
small. The corresponding correction in the underdamped limit is proportional to X 2 Ei ) /(kBT) and 
is not necessarily small. We thus suggest that the effects due to the velocity-dependent friction 
may be of considerable importance in determining the rate of escape of an under- and moderately 
damped Brownian particle from a deep potential well, while they are of minor importance for an 
overdamped particle. 
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I. INTRODUCTION 



Many chemical reactions in a condensed phase can be modeled as an escape of a Brownian 
particle from a potential well. The conventional means for describing this kind of problems 
are the Fokker-Planck equations (FPEs), or master equations, or the (generalized) Langevin 
equations for the reaction coordinate and the conjugated momentum of the particle in the 
external potential. Since the seminal work of Kramers, his approach has further been refined, 
and the influences of different dissipation (collision) mechanisms,— "2*2^ non-Poissonian col- 
lision statistics,-^ non-Markovian effects,- 1 ^ 1 ^ 1 ^ and position-dependent friction^L^L£ on 
the reaction rates have been investigated (see also reviews^!^^) . 

In the present paper, we study the effects due to the velocity dependence of friction — 
Our motivation is as follows. There have been collected evidences in the literature that 
relaxations of linear and angular velocities and kinetic energies in liquids under (highly) non- 
equilibrium conditions are reproduced neither by conventional FPEs with constant frictions 
nor by master equations with constant collision frequencies,— i22i22j212§ so that one has to take 
into account (linear and/or angular) velocity dependence of friction.— i22a2&24i2lii2§i2L2& Indeed, 
if we consider relaxation under equilibrium conditions, then all the relevant energies are of 
the order of k B T {ks being the Boltzmann constant and T being the temperature of the 
bath). Putting aside non-Markovian effects, we can then always introduce a certain effective 
velocity-independent friction. This means that the effects due to the velocity dependence 
of friction are normally small and can safely be neglected. If we consider relaxation under 
non-equilibrium conditions, when initial energies of the particles are much higher than the 
typical bath energies ~ k B T, then the velocity dependence of friction cannot be ignored, in 
general. 

It is not unreasonable to expect that similar effects may become relevant in the evaluation 
of the activated escape rates for underdamped or moderately damped Brownian particles. 
If an underdamped particle is trapped in a deep potential well (Ej, ^> k B T, see Fig. 1) then 
friction is the rate determining parameter: fluctuations must supply the particle with the 
energy ~ E b ^> k B T in order to surmount the barrier. Thus, we have to adequately describe 
the Brownian particle relaxation not only within the thermal energy interval ~ k B T but also 
for much higher energies ~ Eb. In such a case the constant friction approximation might be 
an oversimplification, so that the velocity dependence of friction may become crucial. On 
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the contrary, the effect of the velocity-dependent friction on the the overdamped Brownian 
particle is expected to be minor, since dissipation manifests itself through the diffusion 
coefficient which equals to the integral relaxation time of the equilibrium velocity correlation 
function, which is determined by the typical energies ~ k B T. The present paper is aimed 
at proving and elucidating the above qualitative expectations. 

II. ESCAPE RATES IN RAYLEIGH MODEL 

All the models studied so far within the framework of the theory of activated rate pro- 
cesses treat dissipation either via the FPE with velocity-independent friction or or via mas- 
ter equations with velocity-independent collision frequency-^^ 1 ^^ In order to take into 
account velocity-sensitive dissipation, we must start either from the FPE with velocity- 
dependent friction 28 or from master equations with velocity-dependent collision rates.— 
To avoid introducing phenomenological dissipation models, we use the Rayleigh model of the 
nonlinear Brownian motion, which describes a one-dimensional (ID) relaxation of a heavy 
particle of mass M in the bath of light particles of mass m <C M. Within this model, colli- 
sional relaxation of the heavy particle is described by the master equation with the velocity- 
dependent collision frequency.— 1 ^ The model possesses a small parameter A 2 = m/M <C 1. 
Then, according to van Kampen 2 ^, we can run a systematic expansion of this master equa- 
tion, which allows us to construct the dissipation operator to any desirable order in A. Up 
to the terms of the order of 0(A 4 ), our starting equation reads'^ 2 - 



Here p(x,v,t) is the probability density function in the position (x) - velocity (v) phase 
space, U(x) is an external potential which is schematically depicted in Fig. 1, 



d t p(x, v, t) = {-vd x + M- l U'{x)d v + f (Li + A 2 L 2 ) } p(x, v, t). 



(1) 



87/ 



AV- 1 / 2 



(2) 




is the friction, rj is the number of the bath particles per unit length, 



<t = /3M, (3 



1 



(3) 



k B T 



The dissipation operator in Eq. (TTJ consists of two parts. 

Li = d v v + -<9 2 



a 



(4) 
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is the standard Fokker-Planck operator, while the next-order correction to it is explicitly 
written as 

l 2 = -o v v + a -oy - 2 -di + hy + ±%v + -i^ 4 . (5) 

o a I 6a oa z 

After doing some algebra, we can rewrite Eq. (pQ) in the following compact form: 

dtp(x,v,t) = {-vd x + M- 1 U'(x)d v + ^d v p B (v)S(v)d vP s 1 (v)j p(x,v,t). (6) 

Here 

S(v) = l + \ 2 5 + X 2 p(d v + (q~ <t)v) (d v - qv) (7) 
is the core of the dissipation operator, and 



Pb 



(v) = vV(27r) exp{-™ 2 /2} (8) 



is the equilibrium Boltzmann distribution, and 

P=^ 1 = a — - A > 6 = ± TZ- ( 9 ) 

oa z 4 da 

The standard Kramers FPE is recovered in the limit A 2 — > 0, S(v) — > 1. Eq. ([6]) can be 
solved numerically for any values of parameters (Q. To get a better understanding of the 
effect of velocity- dependent friction on the escape rates, we analyze the problem analytically. 

Applying the projector operator technique, we can reduce 2D Eq. ([6]) to the corresponding 
ID equations in the limit of the high (£ 3> Q) and weak (£ Q) friction, Q being a 
characteristic well frequency. In the former case, the Brownian particle's velocity is a fast 
variable. We can thus introduce the projection operators 

P = Pb(v) J dv, Q = l-P, 

and arrive at the Smoluchowski diffusion equation for the position-dependent probability 
density: 

d tP ( Xl t) = Dd x {(3U'(x) + d x }p(x,t). (10) 

Here the diffusion coefficient, D, equals to the integral relaxation time of the equilibrium 
velocity correlation function in the overdamped limit, J* °° dt (vv(t)). Within the present 
model this correlation function has been evaluated by Hynes, — so that 



and the Kramers rate (in the limit of low temperature and/or high barrier height, /3Eb 3> 1) 
reads 

k = eM~PE b } + 0(A 4 ). (12) 

The frequencies ujq and u b determine the shape of the potential linearized in the vicinity of 
the bottom of the well and at the barrier, 

U(x) » uj {x - x ) 2 /2 and U(x) w E b - u b (x - x b ) 2 /2, (13) 

respectively, see Fig. 1. As expected, the velocity-dependent friction increases the diffusion 
coefficient,— 1 ^ but the effect is of the order of A 2 and is therefore small. 

The underdamped limit (£ <C Q) is much more interesting. We can either switch to the 
action (/) - angle (0) variables and introduce the projectors 

1 f 2n 

P = — dcj>, Q = l-P, 
^ Jo 

or change to the energy variable, E. In any case, we can start from Eq. ([6]) and arrive at 
the energy FPE 



d t p(E,t) = Cd E b(E) Peq (E)d EP ; q 1 (E)p(E,t). (14) 
The core of the energy FPE operator ([6]) is explicitly written as 

D(E) = D Q (E) + D 1 {E)d E + D 2 (E)d 2 E ; (15) 

D (E) = jj{bi(E) + X 2 ((3b 3 (E)/3 - h(E))}, (16) 

D 1 (E) = ^{3b 1 (E)/(3-2b 3 (E)}, (17) 
i a\2 

D,{E) = W h(E). (18) 

The equilibrium distribution reads 

Peq (E) = Zz 1 a_ 1 (E)ex P {-f3E}, Z E = J dEa.^E) ex V {-(3E}, (19) 
and we have introduced the quantities 

a n (E) = [ dx(E - U(x)) n/2 , b n (E) = a n (E) / a_i(£) (20) 
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(integration is presumed over all x for which E > U(x)). The standard Kramers energy- 
diffusion FPEii^>2i is recovered in the limit A 2 —>■ 0. Eq. (|14p can be used to run perturba- 
tive series in A 2 to calculate the mean first passage time, t(E), the inverse of which yields 
the relaxation rate, k (see Appendix). If we again assume that [3E b 3> 1, we get then 

+ 0(A 4 ). (21) 

This is the main result of the present paper. If we neglect the ~ A 2 correction, then we 
recover the Kramers expression for the rate in the underdamped limitr^ 2 - 1 ^ 

k Kr = 2£(3exv{-(3E b }^-. (22) 

On the other hand, the expression in the square brackets in Eq. (12 lj) can easily be evalu- 
ated. For the harmonic oscillator (Fl3l) . for example, it equals 1 + \ 2 j3E b /4:. For the Morse 
oscillator it yields 1 + X 2 /3E b /Q. Therefore, the actual small parameter of the problem in 
the underdamped limit is \ 2 f3E b rather than A 2 . Since Eq. (12T1) has been derived in the 
limit (3 Eh 1, the product \ 2 f3E b can be ~ 1 or higher even for A 2 < 1. In this case, 
evidently, the perturbative expansion of the dissipative operator in A 2 may break down, and 
the description within the Rayleigh master equation 27,29 might be necessary. 

We did not consider the influence of non-Markovian effects on the escape rates. These 
effects normally reduce molecular friction and therefore decrease the escape rates in the 
underdamped limit. Their influence is just the opposite to that of the velocity-dependent 
friction. However, as has been pointed out in,— ^ the non-Markovian effects are of minor 
importance in the case of a deep well, /3Eb 3> 1. It is in this latter case we predict the 
velocity dependence of friction can give rise to a dramatic increase of the escape rate, ~ 
X 2 (E b (3) 1 , despite its influence on the velocity relaxation under equilibrium conditions is 
small, ~ X 2 (E b (3)° < 1. 

III. CONCLUSIONS 

The main result of our paper can be formulated as follows: the effects due to the velocity- 
dependent friction 2 ^ may be of considerable importance in determining the rate of escape 
of an under- and moderately damped Brownian particle from a deep potential well. Similar 
effects due to the (possible) energy dependence of vibrational relaxation rates might also be 



l + A 



2 pa 3 (E b ) 
2>a x (E b ) 
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significant in describing unimolecular reactions under non-equilibrium conditions, like those 
investigated in.— ^ 
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APPENDIX A: CALCULATION OF THE FIRST PASSAGE TIME IN THE UN- 
DERDAMPED LIMIT 

The first passage time in the underdamped limit, t(E), can be evaluated through the 
equation 34 

- 1 = p; q 1 (E)d E D(E)p eq (E)d E r(E) (Al) 

obeying the boundary conditions r'(0) = 0, T{E b ) = 0. The adjoined energy FPE operator, 
p- q l {E)d E D{E)p eq {E)d E , is explicitly defined via Eqs. (|I5|)-(|20||i2fc In order to solve Eq. flAlj) 
perturbatively, note that D (E) possesses contributions both ~ 1 and ~ A 2 , while Di(E) 
and D2(E) are both ~ A 2 . Thus we get 

r{E) = J E dE'^^{l-Y(E')} dE" Peq (E") + 0(X% (A2) 

Y(E) = Peq (E) {D x {E)d E + D 2 (E)d 2 E } P -^. (A3) 

We further assume that the temperature is low, /3Eb ^> 1. Then, taking into account 
the explicit form of D\(E) and D2(E), we see that the leading contributions stemming from 
Y(E), which are of the order of A 2 (_Ef,/3)°, cancel each other. Thus, Y{E) yields contributions 
~ A 2 (-Eb/3) _1 , which can be neglected as compared with those ~ \ 2 (E b j3)° and ~ \ 2 {E h (3) 1 . 
Then, omitting Y(E) in Eq. fl A2[) and evaluating integrals to the leading order, we arrive 
at the expression 

r(E) = ^a-i(^) e W {(3E b }^P- + 0(A 4 ). (A4) 
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Finally, retaining only the leading contribution ~ X 2 (Ei ) f3) 1 into Do(Eb) (|T6l) and inverting 
Eq. (JA4|) one gets Eq. ([HD- 
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FIG. 1: Sketch of the potential energy surface. 
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